################################################################################ 
#
# Counterspeech encouraging users to adopt the perspective of minority 
# groups reduces hate speech and its amplification on social media
#
# Tab D6 - Summary statistics
#
################################################################################ 


rm(list = ls())
setsave = T

################################################################################ 
#  LIBRARIES
################################################################################ 

library(dplyr)
library(Hmisc)
library(readr)
library(tidyr)
library(psych)
library(xtable) 

################################################################################ 
#   DATA AND FOLDER
################################################################################ 

wd = ''
wd_res = paste0(wd, '/results')
wd_data =  paste0(wd, '/data')

data = read.csv(paste0(wd_data, '/dataset_analysis.csv'))
data = data[!(data$user_deleted == "yes"),]


################################################################################ 
#   VARIABLES
################################################################################

outvars = c('orig_tweet_deleted', 'total_tweets_post',
            'tox_dummy_p_key_num', 'tox_dummy_p_key',
            'n_favorites_to_orig_tweet', 'n_retweets_to_orig_tweet',
            'reply_count', 'n_favorites_to_treat_tweet', 'n_retweets_to_treat_tweet')


outvars_name = c('Xenophobic Tweet Deleted', 'No. of Total Tweets',
                 'No. of Xenophobic Tweets', 'Xenophobic Tweet Share',
                 'No. of likes received by xenophobic tweet', 'No. of retweets of xenophobic tweet',
                 "No. of replies to xenophobic tweet", 'No. of likes received by counterspeech tweet',
                 'No. of retweets of counterspeech tweet')

################################################################################ 
#   post treatment outcomes
################################################################################ 

container = data.frame()

for (j in 1:length(outvars)){
  var = outvars[j]
  print(var)
  sum = unname(describe(data[var])[c(2, 3, 4, 5, 8, 9)])
  container = rbind(container, c(outvars_name[j], sum))
}

colnames(container) = c('', 'n', 'mean', 'sd', 'median', 'min', 'max')
# print(xtable(container, digits = c(0,0, 0, 3,3,3,3,3)) ,include.rownames =FALSE)

if (setsave==T){write.csv(container, paste0(wd_res, "/tables/decriptive_stats.csv"))}



